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Abstract 

Wetlands, and especially their littoral zones, are considered to be CH4 emissions hotspots. The recent creation of reservoirs 
has caused a rapid increase in the area of the world's littoral zones. To investigate the effects of water depth and water level 
fluctuation on CH4 fluxes, and how these are coupled with vegetation and nutrients, we used static closed chamber and gas 
chromatography techniques to measure CH4 fluxes in the littoral zone of a large reservoir near Beijing, China, from 
November 201 1 to October 2012. We found that CH4 flux decreased significantly along a transect from open water to dry 
land, from 3.1 mg m~^ h~^ at the deep water site to approximately 1 .3 mg m~^ h~^ at the shallow water site, and less than 
0.01 mg m""^ h~^ in the non-flooded area. Water level influenced CH4 flux by affecting soil properties including soil redox 
potential, soil carbon and nitrogen, and bulk density. The largest emission of all was from the seasonally flooded site after a 
flooding event (up to 21.1 mg m~^ h~^), which may have been caused by vegetation decomposition. Submerged sites had 
greater emissions, while the driest site had lower emissions. Immediately after the monthly measurements had been made, 
we removed the aboveground vegetation to enable an assessment of the gas transportation per unit of biomass. Removal 
of biomass decreased emissions by up to 53%. These results indicated the dominant effect of water depth on CH4 flux 
through effects of soil conditions, plant species composition and distribution. This study suggests that temporally flooded 
wetlands, including littoral zones, contribute significantly to the global CH4 burden. However, the current challenge is to 
capture their spatial extent and temporal variation in the fluxes. 
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Introduction 

Methane is the second most important contributor to total 
greenhouse gas emissions, with a Global Warming Potential of 25 
over a 100-year time-span [1]. Wetlands are considered to be the 
most important CH4 source, thought to account for 24-39% of the 
total global emission of CH4, albeit with a large degree of 
uncertainty [2]. Reservoirs are an important type of wetland, 
whose combined area across the globe has increased in recent 
years to now occupy approximately 5x10^ km^, which is 1/3 the 
size of all natural lakes [3] . One reason for the large development 
of reservoirs is their potential as a clean energy source, although 
uncertainties in their methane emissions have led some scientists to 
question whether they are as clean as people believe [4]. Several 
authors have insisted that the subject of greenhouse gas emissions 
from reservoirs requires more scientific study [5-7] . The Ministry 
of Water Resources of the People's Republic of China has stated 
that China has over 80,000 reservoirs, large and small [8]. 
However, work in China on their greenhouse gas emissions, and 
the extent to which these reservoirs contribute to total emissions, is 
limited [9]. 

The littoral zone of reservoirs features fast and complicated 
material cycles and is reported to be a hotspot for CH4 emissions 
[10,11]. Furthermore, there are indications from elsewhere about 



the importance of the littoral zone: for example, in boreal systems, 
where about 70% of total CH4 emissions comes from the littoral 
zone, even though its area is no more than 24% of the total 
wetland area [12]. Water depth and its diurnal and seasonal 
fluctuations are thought to be the main drivers determining the 
characteristics of the littoral habitat, including plant species 
composition and distribution, soil conditions, and CH4 flux 
[11,13-16]. However, the interaction between water level 
fluctuation, CH4 production, transportation, and emissions is still 
not fully understood [1 7, 18], especially the effect of dynamic water 
fluctuation coupling with vegetation amd nutrients [19,20]. The 
rate of change of the water table, as distinguished from just 
differences in the water table, is known to be a major 
environmental factor controlling the CH4 flux [21]. 

In summary, studies of the CH4 flux from the littoral zone are 
essential for evaluating the atmospheric and climatic impacts of 
reservoirs as energy sources, as well as for a better understanding 
of the biogeochemical mechanism of CH4 flux. Research on the 
mechanisms of how the water level and its fluctuations influence 
CH4 flux is an important foundation for related work aimed at 
modeling and controlling carbon loss through ecosystem manage- 
ment. In this context, we report the seasonal and monthly 
variation of CH4 flux of the littoral zone of Miyun Reservoir, 
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Beijing, China. We hypothesized that water depth and water level 
fluctuation would be the main factors driving the CH4 flux 
variability, together with soil and vegetation. 

Materials and Methods 

This study was authorized by Beijing North Miyun Reservoir 
Eco-agriculture Co. Ltd. It did not harm any protected species, 
and thus no ethics committee was required to authorize the work. 
All data used in this paper were collected in the field by the 
research team, and the raw data are available upon request fi:'om 
the senior author. No data were downloaded firom publicly 
available resources, and thus no permission was required. 

Study area 

The research was carried out at Miyun Reservoir (40°29'N, 
116°50'E), which is located in the northern mountainous area 
near Beijing, China. It was built in 1960 with a maximum water 
area of 188 km^. Its catchment is characterized by warm 
temperate semi-humid monsoonal climate with an annual average 
air temperature of 10.5°C, maximum air temperature of 38°C, 
and a minimum of — 18°C. The reservoir is normally covered by 
ice from the middle of November to the end of March. The 
growing season is from April to November. The annual average 
precipitation is close to 600 mm, of which 80% is concentrated 
from July to August [22]. The annual change in the water level is 
1-5 m because of rainfall and water supply for domestic use. The 
area between the highest and lowest water level from 1984 to 2005 
was 84 km^ [23]. In the summer of 2012, when the work was 
carried out, continuous heavy rain in July caused a sudden water 
level increase of one meter, and part of the littoral vegetation was 
inundated. 

We divided the littoral zone into five areas based on water level 
(Figure 1). Sites were selected ranging from locations in open water 
to the dry area on higher ground, to provide five contrasting 
environments: (i) deep water area (DW); (ii) shallow water area 
(SW); (iii) seasonal (August and September) flooded area (SF); (iv) 
'seasonally flooded control' (SFC) area, which was 500 m away 
from SF, had the same plant species and similar soil carbon/ 
nitrogen content as SF, but escaped the flood in August and 
September because of its 1-m-higher elevation; and (v) permanent 
non-flooded area (NF). Details of the water depths in each of these 
areas are shown in Figure 2B. Three typical plant communities in 
each area were selected, and Table 1 shows the dominant species 
in different seasons in all zones. For more details on biomass and 
soil, see Figures 2 and 3. 

CH4 flux measurements 

Methane flux was measured in November 201 1, then May, July, 
August, September and October 2012. The experiment of three 
plots at site SFC was carried out just after the flooding and during 
the time when the water level dropped from August to October 
2012. In order to reduce uncertainty in the average daily flux, 
sampling to capture any diurnal variation was performed at three- 
hourly intervals (local time: 6, 9, 12, 15, 18, 21 and 24 h). Each 
plot had four replicates located within three meters from each 
other. To eliminate disturbance to the soil, wooden access 
platforms were built. 

The static opaque chamber technique was used to determine 
the CH4 flux [24]. The chambers were made of stainless steel 
(volume: 125 L; surface area: 0.25 m^) and covered with 
polyethylene foam to avoid any warming effect inside the 
chamber. An internal chamber (volume: 200 L; surface area: 
0.25 m^) could be added if plants were tall. Two fans were built 




Figure 1. Experimental design. WL: water level. The sites are 
grouped at different heights. DW: deep water site; SW: shallow water 
site; SF: seasonally flooded site; SFC: 'control site' for the seasonally 
flooded site; NF: non-flooded site. A, B and C denote samples from 
different vegetation types within each height band (see Table 1 for 
information on plant species). There were four replicates in each case, 
repeatedly sampled six times (also repeatedly sampled seven times in a 
day) in the year. For more details on water depth and other 
environmental parameters, see Figures 2 and 3. 
doi:10.1371/journal.pone.0094275.g001 

into the chamber for air mixing. Four gas samples (200 mL each) 
were taken using 100-mL polypropylene syringes at 15-min 
intervals over a 45-min period after enclosure, and stored in 
500-mL plastic and aluminum membrane gas sampling bags. The 
concentration of CH4 was analyzed within one week by gas 
chromatography (7 890 A, Agilent, California, USA) equipped with 
a flame ionization detector (FID). Gases were separated with a 
column (3 m, 3.2 mm) packed with Porpak (80/ 100 mesh). The 
temperatures of the oven, injector, and detector were 70°C, 20°C, 
and 200°C, respectively. The flow rate of the carrier gas (N2) was 
25 mL min~\ Standard CH4 gas (2.03 ppm in air, China 
National Research Center for Certified Reference Materials, 
China) was used for precision verification for CH4 concentrations. 
The coefficient of variation was below 1%. The flux of CH4 was 
calculated following Chen et al. [25]. In order to determine the 
effect of aboveground vegetation on ffux, one more ffux 
measurement was taken at 9 am the following day (after seven 
times sampling for diurnal variation), with aboveground plant 
material removed. Chambers were reset into new positions near 
the old positions each sampling month. All positions at each site 
were within an area of 20 m^, but not so close to each other to 
cause artifacts in the data through (for example) changes in the 
local hydrology. 

Environmental factors 

In order to analyze the effects of environmental parameters on 
fluxes, the following factors were taken into account: water level, 
dissolved oxygen (DO) in water, soil total carbon (TC) and 
nitrogen (TN), soil bulk density, biomass, and air temperature. 

Water level was measured after gas sampling at DW, SW and 
SF (when SF had standing water in August and September 2012). 
At site SF (when there was no standing water in November 201 1, 
May, July and October 2012) and SFC, a 1-m PVC tube was 
inserted vertically into the soil under the chamber after all monthly 
gas sampling was complete, allowing two hours for the water level 
to equilibrate before measuring the level. The water table of site 
NF was calculated according to the elevation measured by a 
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Figure 2. Monthly variation of air temperature, water level and 
biomass of each site. Days between dotted lines was the high water 
level period and thus the seasonally flooded site (SF) was under water. 
DW: deep water site; SW: shallow water site; SF: seasonally flooded site; 
SFC: 'control site' for the seasonally flooded site; NF: non-flooded site. A 
negative value of water level indicates the ground water depth. Error 
bars represent SE. 

doi:1 0.1 371 /journal. pone.0094275.g002 
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Figure 3. Soil properties (mean ± SE) of each site. Bars with 
different letters indicate a significant difference at p<0.05. TC: total 
carbon, TN: total nitrogen. DW: deep water site; SW: shallow water site; 
SF: seasonally flooded site; SFC: 'control site' for the seasonally flooded 
site; NF: non-flooded site. 
doi:10.1371/journal.pone.0094275.g003 

Global Navigation Satellite System receiver (BLH-L90, Daheng 
International, China). DO in water was measured during the gas 
sampling by a handheld multi-parameter meter (Professional Plus, 
YSI, USA), after flooding. 

Soil samples at site DW, SW, SF and NF were collected from 
three different layers (0-10 cm, 10-20 cm, and 20-30 cm below 
ground) at each replicate location in November 2011, except site 
SFC in October 2012. After air-drying and grinding (passing 
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Table 1. Dominant plant species at each plot in different months. 





Site 




Nov 201 1 


May 2012 Jul 2012 


Aug 2012 Sep 2012 


Oct 2012 


DW 


A 


Echinochloa olonum^^ 


Myrioptiyllum sp. 


Trapa'^^ sp. 






B 


no vegetation 










C 


Typha angustifolia'^^ 








SW 


A 


Xanthium sibiricum^ 


Scirpus planiculmis^^ 


Echinochloa colonum^^ 






B 


Setaria viridis^ 


Bidens pilosa^ 


Echinochloa colonum^^ 






C 


Zeo mays^ 


Polygonum la pattii folium^ 




Typha angustifolia^^ 


SF 


A 


Xanthium sibiricum 


Cirsium setosum 


Cirsium setosum^ 


Cirsium setosum 




B 


Setaria viridis 


Hemarthria altissima 


Hemarthria altissima^ 


Hemarthria altissima 




C 


Zea mays 


Polygonum lapathifolium 


Polygonum lapathifolium^ 


Polygonum lapathifolium 


SFC 


A 


# 




Cirsium setosum 






B 


# 




Hemarthria altissima 






C 


# 




Zea mays 




NF 


A 


Xantliium sibiricum 










B 


Setaria viridis 


Artemisia argyi 








C 


Zea mays 









# indicates no data; 

DW: deep water site, SW: shallow water site, SF: seasonally flooded site, SFC: 'control site' for seasonally flooded site, NF: non-flooded site; 
A, B, C indicates sample plot with different vegetation; 

Species with aerenchyma are denoted ^ species that are emergent are denoted ^ 
doi:1 0.1 371/journal.pone.0094275.t001 



through a 100 mesh sieve), TC and TN were analyzed using an 
elemental analyzer (vario MACRO cube, Elementar, Germany). 
Soil bulk density was measured following Chinese national 
standards NY/T 1121.4-2006 [26]. 

To determine which species had aerenchyma tissues, thin hand 
sections were made of roots and stems (photomicrographs were 
made). The aboveground biomass of every replicate in the 
chamber was weighed after drying at 80°C to constant mass. 

Diurnal air temperature was measured by a digital thermometer 
(JM624, Jinming, China) at the start and end of each gas sampling 
at every plot. 

Statistical analysis 

Reported daily fluxes were obtained by averaging the means of 
the three-hourly values. To estimate how much CH4 was 
transported by vegetation and to calculate the differences between 
species, plant-mediated flux was estimated on the basis of the 
difference between the flux before and after vegetation removal. 
We followed Kankaala et al. [27] to derive a mass-based plant 
transport index, PTI: 



flux with aboveground vegetation —flux without aboveground vegetation 
biomass 

The Kruskal-Wallis H test was used to test for spatial difference 
of flux and the Mann-Whitney U test for further multiple 
comparisons. Spatial differences of soil properties were analyzed 
with one-way ANOVA, and then LSD for multiple comparisons. 
A logio transformation was used to show the correlation between 
water depth and CH4 flux. A value of 0.5 was added to the flux 
data before transformation to make sure that all of the data were 
positive. Spearman analyses were used for correlations between 



flux and environmental factors. All the analyses above were 
performed using IBM SPSS Statistics (19.0, IBM, USA). Charts 
were made using SigmaPlot (11.0, SYSTAT, USA). 

Results 

Spatial and seasonal variation of CH4 flux 

CH4 fluxes at the five sites showed significant differences 
(n = 324, p<0.05; Figure 4). CH4 fluxes at permanently flooded 
sites, both the deep and shallow water sites (DW: 3.1 ±0.5 mg 
m h"^; SW: 1.3±0.2 mg m h ), and the seasonally flooded 
site (2.1 ±0.4 mg m ^ h ^) were significantly higher than non- 
flooded sites. The seasonally flooded site SF, which was flooded for 
two months, emitted higher (1.6 times) CH4 than the permanent 
flooded site SW with shallow water. The seasonally flooded site SF 
was dramatically higher than its control site, SFC (the flux was 247 
times higher). Although the never-flooded site NF presented as a 
weak sink (—2.7 ±9.2 10 ^ mg m ^ h ^), while site SFC was a 
weak source (8.7 ±8.9 10~^ mg m~^ h~^), there was no significant 
difference between them (n=144, p>0.05). The average CH4 
emission from all the flooded sites (DW, SW and SF) was about 
1 120 times higher than the non-flooded sites (SFC and NF). 

There were different seasonal patterns among different sites 
(Figure 5). CH4 flux at site DW and SW continuously increased 
from November 2011 to August 2012, and thereafter decreased 
gradually. CH4 flux at site SF increased slightly from November 
2011 to July 2012. After flooding it rose sharply (from 0.05 mg 
m ^ h ^ to 6.4 mg m ^ h ^) to a high peak, and then remained 
the highest emitter among all sites until the autumn. The seasonal 
pattern of site NF was totally contrary to sites DW, SW and SF; 
although the amplitude of flux was narrow. It decreased from 
autumn to the following summer, and then increased gently. CH4 
flux at site SFC was as low as site SF before it was flooded. 
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DW SW SF SFC NF 



Site 

Figure 4. CH4 flux (mean ± SE) at different sites. Bars with 
different letters indicate a significant difference at p<0.05. Positive 
values of flux indicate a CH4 source. Inset in the top right-hand corner 
shows SFC and NF on a different scale. DW: deep water site; SW: shallow 
water site; SF: seasonally flooded site; SFC: 'control site' for the 
seasonally flooded site; NF: non-flooded site. 
doi:1 0.1 371 /journal. pone.0094275.g004 




-2 1 1 1 1 1 1 1 1 1 — ' — I 1 ' I ' 

NDJFMAMJJASO 
2011 2012 

Figure 5. Monthly variation of CH4 flux (mean ± SE) at different 
sites. A positive value of flux indicates a CH4 source. Inset in the top 
right-hand corner shows SFC and NF on a different scale. Days between 
dotted lines was the high water level period and thus the seasonally 
flooded site (SF) was under water. DW: deep water site; SW: shallow 
water site; SF: seasonally flooded site; SFC: 'control site' for the 
seasonally flooded site; NF: non-flooded site. 
doi:10.1371/journal.pone.0094275.g005 



Effects of environmental factors 

Considering all the cases where the water depth was 2 m or less, 
the water depth was positively correlated with CH4 flux (n = 324, 
p<0.01; Table 2). When there was standing water (water depth 
>0 cm; Figure 6), the habitat always emitted more CH4 than 
when there was no standing water (water depth <0 cm). 

Soil properties varied according to the zone: TC and TN were 
highest at site DW and lowest at site NF, while bulk density 
showed an opposite trend (n = 60, p<0.05; Figure 3). Soil TC and 
TN were positively correlated with flux at all sites (n = 324, 
p<0.01; Table 2), as well as individually at sites DW and SFC 
(n = 72, p<0.01). Bulk density was negatively correlated with flux 
at the whole-site scale (n = 324, p<0.01; Table 2), and no 
significant correlations were found at separate sites (n = 72, 
p>0.05). 

There were large differences in CH4 flux before and after the 
removal of aboveground vegetation at all sites (Figure 7). After 
vegetation removal, fluxes at sites DW and SW decreased by 
approximately 50% (47% and 53%, respectively), while the 
contribution of transportation at site SF was much smaller (6%). 
For the two emergent sites, SFC and NF, where fluxes were very 
small, opposite patterns were found. After the plants were 
removed, emissions at site SFC increased, while at site NF they 
decreased. The plant transportation index differed among the five 
sites and 14 plots (Figure 8). It showed seasonal variation, with the 
highest transportation efficiency found in spring or summer at 
most of the sites. 

CH4 emissions showed different patterns with different plant 
species after flooding at site SF (Figure 9). High CH4 emissions 
(from 9.6 mg m ^ h ^ to 21.1 mg m ^ h ^; average: 
14.2 ±2.5 mg m ^ h ^) were recorded at the plot with Cirsium 
setosum after the flooding (August; flooded for 20 days), which was 
6.8 and 5.1 times higher than the site with Polygonum lapathifolium 
(from 1.6 mg m ^ h ^ to 2.6 mg m ^ h ^; average: 2.1 ±0.2 mg 
m ^ h ^) and Hemarthria altissima (from 2 mg m ^ h ^ to 4. 1 mg 
m ^ h ^; average: 2.8 ±0.5 mg m ^ h ^), respectively. At the 
same time, the lowest water DO (30%-33%) was observed at the 
plot with C. setosum when the flux was the highest. Water DO was 
negatively correlated with flux at all flooded sites (DW: n = 72, 
p<0.01; SW: n = 72, p<0.01; SF: n = 24, p<0.01; Table 2). 



The vegetation biomass began to increase in July and peaked in 
September to October (Figure 2C). There was a significantly 
negative correlation (n = 324, p<0.01; Table 2) between flux and 
biomass when all sites were taken as a whole. However 
significantly positive correlations were observed at sites DW and 
SF (n = 72, p<0.01; n = 72, p<0.05). No significant correlations 
were found at the other three sites (n = 72, p>0.05). 

Air temperature positively correlated with CH4 flux at sites DW 
and SW (n = 72, p<0.01; Table 2), but was negatively correlated 
at site NF (n = 72, p<0.05). No significant correlations were 
observed at sites SF and SFC (n = 72, p>0.05; n = 36, p>0.05). 

Discussion 

High CH4 emissions from the littoral zone 

The average CH4 emission rate (1.3 mg m ^ h ^) from the 
littoral zone of Miyun Reservoir was 6.5 times higher than the flux 
from the open water area (0.2 mg m ^ h ^) [28]. The low 
emissions from open water are similar to those from reservoirs at 
higher or lower latitudes, but in those cases the littoral zone has 
not been investigated [29,30]. Reservoirs differ enormously in the 
percentage of their area that can be deemed 'littoral'. In the case 
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Figure 6. Relationship between water depth and CH4 flux. 
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Table 2. Correlation coefficient R between CH4 flux and environmental parameters. 



Site 


Water level 


Soil TC 


Soil TN 


Soil bulk density 


Biomass 


Water DO 


Air temperature 


all sites 


0.75^^(324) 


0.62^^(324) 


0.56^^(324) 


-0.53^^(324) 


-0.26^^(324) 


-0.28^^(168) 


0.25^^(324) 


DW 


0.38^^(72) 


0.45^^(72) 


0.43^^(72) 


-0.02(72) 


0.42**(72) 


-0.31^^(72) 


0.63**(72) 


SW 


0.42^^(72) 


-0.06(72) 


-0.05(72) 


-0.03(72) 


0.08(72) 


-0.32^^(72) 


0.65^^(72) 


SF 


0.74^^(72) 


0.19(72) 


0.12(72) 


0.07(72) 


0.24^(72) 


-0.5r(24) 


0.08(72) 


SFC 


0.47^^(36) 


0.57^^(36) 


0.57^^(36) 


0.24(36) 


-0.31(36) 


# 


0.10(36) 


NF 


-0.02(72) 


-0.14(72) 


-0.15(72) 


0.05(72) 


-0.03(72) 


# 


-0.28^(72) 



# indicates no data; 

Numbers in () indicate the sample size; 

TC: total carbon, TN: total nitrogen, DO: dissolved oxygen; 

** indicates significant correlation (P <0.01), ^ indicates significant correlation (P<0.05). 
doi:1 0.1 371/journal.pone.0094275.t002 



of Miyun, in the period from 1984 to 2005, the average area of the 
water surface was 103±15 km^, while its annul difference was 
9.5 ±2 km^ [23]. Although the seasonally flooded area is as high as 
8.5% of the total area of the reservoir, it explains how the littoral 
component of emissions is as much as 40% of the total. The high 
level of emissions from the littoral zone emphasizes its status as a 
methane hotspot in this case, but not all reservoirs will have a 
littoral zone comprising such a high percentage of the total. 
Although the present study involving 15 typical habitats, there are 
still uncertainties in the estimation because of high spatiotemporal 
flux variation and environmental heterogeneity, as demonstrated 
by the complex effects of plant, soil and temperature on fluxes 
found in this research. Furthermore, the dynamics of the littoral 
zone exacerbates the uncertainty. Nevertheless, emissions from the 
littoral zone are likely to form a substantial part of the reservoir's 
methane budget, and deserve further research [10,31]. 




-1 1 1 1 1 1 

DW SW SF SFC NF 

Site 

Figure 7. Contribution of aboveground plants to CH4 flux 
(mean ± SE) of each site. A positive value of flux indicate a CH4 
source. A positive rennainder nneans nnethane is transported fronn the 
vegetation to the atmosphere. Inset in the right-hand center shows SFC 
and NF on a different scale. DW: deep water site; SW: shallow water site; 
SF: seasonally flooded site; SFC: 'control site' for the seasonally flooded 
site; NF: non-flooded site. 
doi:1 0.1 371 /journal. pone.0094275.g007 



Water level fluctuation and increased CH4 emissions from 
the littoral zone 

After 20 days of flooding, CH4 emissions from site SF increased 
sharply to the highest level among all sites, and this high rate lasted 
to October (Figure 5). In contrast, a decreasing of CH4 emissions 
was reported at a site in Finland and was explained by diminished 
plant-mediated transportation and decreased carbon supply from 
growing plants to the rhizosphere after plants were totally 
submerged [16]. In another study, different patterns were shown 
when flooding occurred in different seasons; CH4 emissions were 
negligible after flooding in autumn or winter, but increased sharply 
after flooding in summer [32]. Based on the present work, we 
suggest that anaerobic plant decomposition, which apparently 
occurred at site SF after flooding, is the main reason for high 
emissions. Elsewhere, it has been shown to be the cause of CH4 
loss from tropical reservoirs [33]. Our situation was similar, but 
one in which the effect was more rapid due to the herbaceous 
environment. 

Different fluxes among flood-tolerant and -intolerant plants 
were observed. In the case of the flood-intolerant C. setosum, 
flooding caused a sharp increase in emissions, but for the flood- 
tolerant P. lapathifolium and H. altissima, under the same water 
depth, no sharp increase was seen. This phenomenon might be the 
reason for the difference between the aforementioned research in 
Finland and that of our own; their research was carried out at 
places where Carex is the main vegetation, which does not rot 
immediately after flooding [16]. Different responses of plant 
species to flooding help to explain the spatial variation, and this 
could represent a potential approach to reducing methane 
emissions, i.e., through careful vegetation management. 

The lowest water DO found in the plots with flood-intolerant 
species (Figure 9) may be the result of oxygen consumption during 
plant tissue decomposition [34]. Variations in CH4 emissions have 
sometimes been attributed to variations in methane oxidation [35] . 
The observed negative relationship between DO and CH4 flux 
(Table 2) asserts the importance of available DO and plant 
decomposition on CH4 loss. A dramatic rise in CH4 emissions 
appeared when the DO fell to approximately 30%, suggesting that 
30% might be the lower bound of usable DO concentration for 
methanotrophs. However, more work on this aspect is needed, 
since GH4 oxidation is a complex process affected by many factors, 
including microorganisms, temperature, water content, and redox 
potential [36,37]. It is therefore probable that the DO threshold 
will vary across different habitat types. 

It is conceivable that wetlands featuring a littoral zone 
constituting a relatively large proportion of the total area, as well 
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Figure 8. Above-ground plant transport index (mean ± SE). AF (mg CH4 h~^) is the flux with vegetation minus the flux without 
vegetation, i.e. the flux of methane that can be attributed to the presence of the vegetation. B (kg m~^) is biomass. A positive flux means methane is 
transported from the vegetation to the atmosphere. Days between dotted lines was the high water level period and thus the seasonally flooded site 
(SF) was under water. DW: deep water site; SW: shallow water site; SF: seasonally flooded site; SFC: 'control site' for the seasonally flooded site; NF: 
non-flooded site. 

doi:10.1371/journal.pone.0094275.g008 



as frequent water level fluctuations, might produce especially high 
emissions when the majority of the plant community comprises 
flood-intolerant species. It is further conceivable that climatic 
variability, i.e., an increased risk of flood and drought [38], will 
exacerbate methane production because of CH4 emissions caused 
by unusual water level fluctuations, and those associated with re- 
flooding. For example, these conditions occurred in China after an 
extreme drought in 2010-2011, which caused the drying out of 
numerous lakes, giving them the appearance of grasslands. 



exemplified by a well-documented case at the largest freshwater 
lake in China, Poyang, which has an area of 3,150 km^ [39]. 

Effects of water depth on the flux through soil and plants 

Water depth determines soil and plant conditions, which in turn 
influence CH4 flux through different mechanisms. However, we 
found these effects to be sometimes limited by the water depth 
itself 
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Figure 9. Relationship between dissolved oxygen (DO) of water 
and flux of the seasonally flooded site (SF). 
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Flooding decreases the soil redox potential [40-42], which 
affects both methanogenesis activity and gas transfer through plant 
tissues [43]. Besides, the anaerobic conditions inhibit aerobic 
respiration and thus allow more organic matter accumulation 
[44,45], which was the likely reason for the soil nutrients pattern 
found in our research (Figure 3). Carbon and nitrogen compounds 
form the substrate for GH4 production [46,47]. Positive correla- 
tions are usually observed between soil nutrients and CH4 
emissions [48,49]. However, the insignificant correlations at SW, 
SF and NF suggested that the effect of soil nutrients was sometimes 
weak, and even dominated by other factors. A similar result has 
also been shown at a grassland site in Germany, where flux was 
found to be positively correlated with soil moisture but not 
nitrogen fertilizer application from 0 to 450 kg N ha ^ yr ^ [50] . 
High accumulation rates of organic matter forms a soil that has a 
low bulk density [51] and high porosity [52], which might 
accelerate CH4 diffusion in soil pores, leading to methane release 
[53]. Our results agree with the literature, showing that the water 
level influenced CH4 production through the soil redox potential 
and availability of organic substrates. 

Plants are another important reason for the spatial pattern of 
flux. An apparent interspecific difference in the plant-mediated 
flux was found at site DW (Figure 8A), i.e., the flux transported by 
emergent plants was much higher than that of totally submerged 
plants. Higher emissions from emergent plants have also been 
observed at a Qinghai-Tibetan Plateau wetland, although they did 
not reduce the flux through the water surface [54]. Emergent 
macrophytes with a well-developed vascular system and with 
aerenchyma can be an important pathway for CH4 transportation 
from sediment to atmosphere [55]. Although submerged plants 
have also developed aerenchyma, even a thin film of water over 



the leaves and stems might inhibit gas release from plant tissues 
since the diffusive velocity of gases in water is much slower than in 
air, by a factor of 10"^ [56]. In such cases, emissions will be 
suppressed until the time when ebullition begins. Further 
interpretation in this regard came from a study in which it was 
found that CH4 emissions in a sedge-dominated zone decreased 
significantly, when the flood level was high enough to submerge 
the venting structures of the plants [16]. 

Much variation in plant-mediated CH4 emissions among the 
five sites was observed (Figure 7), decreasing by 99% along the 
water depth gradient from water to dry land. Around 50% of the 
CH4 emissions was transported by plant tissues at the permanently 
flooded sites (DW and SW), as was also found in a marsh with 
Spartina altemiflora [57]. Plant-mediated fluxes from the same 
species growing under different water depths were different, e.g., 
sites SF and SFC (Figures 9C and 9D). This shows that, although 
the transportation ability was different among plant species, it was 
nevertheless highly dependent on flooding. 

Besides the spatial pattern, plant development is also considered 
as an explanation for the seasonal variation of CH4 flux. There 
were positive correlations between flux and the biomass of sites 
DW and SF (Table 2). Gas transportation through aerenchyma is 
the most likely explanation [55,57]. Besides, we also calculated 
methane flux per kilogram of biomass as an index of transpor- 
tation efficiency (Figure 8). The arched patterns of plant 
transportation at sites DW and SW suggest that vegetation activity 
might be correlated with plant transportation capacity [58-60]; 
and the high level of emissions in summer were caused not only by 
the high biomass, but were also influenced by the high 
transpiration rate of vegetation (methane dissolved in plant water 
can be transported from soil to the leaves, and hence reach the 
atmosphere). The interaction between biomass and transportation 
efficiency might be a reason for the lack of correlation between 
biomass and emissions at sites SW, SFC and NF (Table 2). Taking 
all sites as a whole, a negative correlation between flux and 
biomass was observed (Table 2), but this relationship is affected by 
the fact that the lowest CH4 emissions and highest biomass were 
found at the driest site. 
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